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of the masses of possible heavy tetraquarks using the color-magnetic interaction with the flavor 
symmetry breaking corrections. 
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I. INTRODUCTION 

o : 

The past several years have witnessed a renaissance of charmonium spectroscopy. A few charmonium-like new 
resonances X(3872) X(3940_) 0, Y(3940) 3] and Y(4260) were discovered experimentally. Because of its very 
narrow width, proximity to DD* threshold and special decay patterns, Swanson proposed that X(3872) is mainly a 
D°D°* molecular state with some admixture of pJ/if) and tuJ/tp components 0,0]. But the possibility of X(3872) as 
a conventional 1 ++ charmonium is still not excluded although its mass is smaller than the quark model prediction 
Q. Li suggested that X(3872) could be a charmonium hybrid meson However, both lattice QCD simulation 
I , and flux tube mode ^(j indicated that the hybrid charmonium lies in the range 4200-4400 MeV. 

Maiani et al. argued that the X(3872) is a bound state of a pair of diquark and anti-diquark ^lj. With the 
\Q ] spin-spin interaction and X(3872) mass as input, they predicted a 2 + cqcq tetraquark around 3952 MeV, which might 
■ be identified to the X(3940) This 2 + tetraquark state can decay to J ftp plus a light vector meson, and DD via 
' D-wave. 

The multi-quark system was studied many years ago using a quark model with color-magnetic (CM) interaction 
[l2l Il3l Ihl Hfij . Recently, Richard and Stancu used the model to estimate masses of qcqc tetraquarks [3,0,0]. In 
the complete qcqc configuration space, they diagonalized the 6x6 Hamiltonian matrix. The CM interaction parameter 
"--^ ' was determined by the charmed baryon masses. They found the mass of the isoscalar state is around 3910 MeV, which 
; lies well above X(3872) and DD* threshold. 

In a previous work on light tetraquarks fl9| , we found that the strength of CM interaction in tetraquarks is quite 
different from that in conventional qq mesons. The strength of CM interaction in charmed baryons is also different 
from that in the tetraquark system because tetraquarks and charmed baryons have different spatial size. Hence the 
extraction of the CM interaction parameter in the tetraquark system from the charmed baryon is not justified. 

In this work, we identify X(3872) as a tetraquark and use its mass as input to fix the CM strength in the cqcq 
system. In Section[n]we present the simple model Hamiltonian with the color-magnetic interaction for the tetraquark 
system and construct the color-spin-flavor wave functions of + , 1 + and 2 + tetraquarks. Section ITTT1 is the numerical 
analysis. The last section is a short summary. 
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II. HEAVY TETRAQUARK MASSES 

For the tetraquark, the model Hamiltonian reads 

H = J2 m i + H cm (1) 

i 

where to$ is the i-th constituent quark mass. Hg iu is the color-spin interaction Hamiltonian which is derived from 
one-gluon exchange in MIT bag model |12l fl3l l20l |21| : 
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Hcm = - / , VijXi ■ Xjai ■ <Tj (2) 
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where Si is the quark spin operator and \ the color operator. For the anti-quark, \q — —A* and Sq = —a*. The 
values of the coefficients Vij depend on the multiquark system and specific models. In general, takes different 
values for qq, qqq and qqqq systems Here we discuss the heavy tetraquarks, the in the heavy quark system 
are also different from the light quark system. We use the following convention 



where v depends on the wave function of the multiquark system. With given flavor context qi<Z29394, the expression 
of the CM matrix element between two states |fc) and |^) reads 

Vcm (91929394) = Via (?i 32) + Vi 3 (qiq 3 ) + Vuiqiqi) + ^23(9293) + ^24(9294) + ^34(9394) , (4) 

where Vij = (k | XiXj<7i<fj | I). The base states and \l) are chosen to be the eigenstates of SU(6) CS . Since the \k) 
and \l) are not eigenstates of the CM interaction, we need calculate every individual term in Eq. Q and diagonalize 
the {k I Vcm I 1} matrix in the base space in order to obtain the physics mass. 

We use the qq ® qq to construct the color-spin wave function of + , 1 + and 2 + heavy tetraquarks. A particular 
multiquark configuration is denoted as [Dq, D 3c , S, N), where Dq and D 3c are SU(6) CS color-spin and SU(3) C color 
representations of the multiquark system respectively. S is the spin of the system and N is the total number of quarks 



and antiquarks. The SU(6) CS eigenstates of the + , 1 + and 2 + heavy tetraquark are listed below: 

|l,l c ,0,4) = 75|21, 6 , 1,2)® |2T,6 c ,l,2) + yi|21,3 c ,0,2)® |2l,3 c ,0,2) (5) 

|405, l c , 0, 4) = Jl\21, 6 C , 1, 2) ® |2l, 6 C , 1, 2) - W~|21, 3 C , 0, 2) ® |21, 3 C , 0, 2) (6) 

|l,l c ,0,4) = 7||15,3 C ,1,2)®|T5, 3,1,2) + a/||15 5 6 c ,0,2> (g |15,6 C , 0,2) (7) 

|189, l c , 0, 4) = y||15, 3 C , 1, 2) ® |T5, 3, 1,2,)- 7||15, 6 C , 0, 2) <gi |15, 6 C , 0, 2) (8) 

|35, l c , 1, 4, ) = \2l,6 c , 1,2,)® I2X6C, 1,2,) (9) 

|35, l c , 1,4) = |15, 3 C , 1, 2, 6/) ® 1 15, 3 C , 1, 2) (10) 

|35, l c , 1, 4) = y||21, 3 C , 0, 2) ® |15, 3 C , 1, 2) - 7||21, 6 C , 1, 2) ® |T5, 6 C , 0, 2) (11) 

|280, l c , 1, 4) = y||21, 3 C , 0, 2) ® |15, 3 C , 1,2,) + ^|21, 6 C , 1,2,)® |T5, 6 C , 0, 2, ) (12) 

|35, l c , 1, 4) = ^115, 3c, 1, 2) ® |2T, 3 C , 0, 2) - ^||15, 6 C , 0, 2) ® \21, 6 C , 1, 2) (13) 

|280, l c , 1, 4) = y||15, 3 C , 1, 2) ® |2T, 3 C , 0, 2) + 7||15, 6 C , 0, 2) ® |2T, 6 C , 1, 2) (14) 

|405, lc, 2, 4) = |21, 6, 1, 2) ® |21, 6, 1, 2) (15) 

|189,l c ,2,4) = 1 15, 3, 1,2)® 1 15, 3, 1,2) (16) 



Here we focus on the + , 1 + and 2 + states without radial and orbital excitations. Their spatial wave functions are 
symmetric. Therefore their color-spin wave functions of flavor-symmetric diquark {qq} or {qq} inside tetraquarks are 
antisymmetric. Those flavor-antisymmetric diquarks [qq] or [qq] have symmetric color-spin wave functions. We do 
not assume that quarks in the tetraquark really form diquarks with strong correlation. Instead we simply use the 
qq ® qq to construct the convenient base wave functions. On the other hand, one can also start from the qq ® qq basis. 
In order to calculate the VcAf (99), we must do some recoupling from the [qq, Dq, D 3c , S, 2) ® [qq, Dq, D 3c , S, 2) basis 
to the [qq, Dq, D 3c , S, 2) ® [qq, Dq, D 3c , S, 2) basis. Details can be found in Ref. [l9|. The recoupling coefficients are 
listed in the Appendix. 

III. NUMERICAL RESULTS 

In order to estimate the heavy tetraquark mass, we need determine the constituent quark mass and the value of the 
parameter v in Eq. J2J for the heavy tetraquark system. Using the experimental values of J/tp, Vc, D*, D, D* s , D s , T 
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masses, we use Eq. to extract the constituent quark mass: 



'M(J/i>) 


= 2m c 






= 1534MeV 


M{ Vc ) = 


2m c - 




m u 


= 437MeV 


M (D*) = 


= m„ + 




m s 


= 542MeV 


M(D) = 


m„ + 




m b 


= 4730MeV 


M(D* S ) = 


= m s + 




Vcu 


= 23.4MeV 


M(D S ) = 


= TOs + 


m c -16M^)(^) 


Vcs 


= 29.3MeV 


k M(T) = 


2rrib + 




Vcc 


= 67.7MeV 



where v ca , v cS and v c£ are merely the CM interaction parameters in Eq. @ for the cc, cu, cs meson system respectively. 
With the assumption that X(3872) is a qcqc tetraquark, the VcAi(qcqc) reads 

V CM {qcqc) = CcV 12 + Via + ( C V 14 + CcV 23 + £V 24 + ( C V 34 (18) 

where ( c = After diagonalizing the full matrix (k \ VcM(qcqc) \ I), we obtain CM interaction eigenvalue and 
mass. The lowest eigenvalue is 

V C m = -15.9uq , 

where vq — v qC qc- In the heavy quark limit Cc = rr 4 — > 0, the CM energy of qcqc system is dominated by the qq 
subsystem VcAi(qcqc) — > V13 = —16vq. The wave function of this lowest-lying tetraquark state will be |cc35, l c , 1, 2) ® 
|ggl, l c , 0, 2) with charge conjugate parity C = — . This state is strongly coupled to J/tp and a light pesudoscalar. 
Therefore it fall parts very easily and has a very broad decay width. Thus it is unlikely to observe it experimentally. 
The second lowest-lying tetraquark state will be |cc35,8 c , 1,2) <X> \qq35,8 c , 1,2) with C = +. Wc identify it to the 
observed X(3872) state: 

Vcm(X(3872)) = -A.lv Q , 

M(X(3872)) = 2m c + 2m u - 4.1vq . 

In the heavy quark limit ( c = ^ — > 0, VcM(qcqc) — > Vi 3 = — 4uq. Using the experimental value of X(3872) mass as 
input, we get 

v Q w 17.1MeV . 

We will use this vq to estimate the mass of other qqQQ and qQqQ tetraquarks. 

The CM strength parameter (v) depends on the constituent quark mass and spatial wave function of quark model. 
Since we consider only the + , 1 + and 2 + heavy tetraquarks without radial and orbital excitations, their spatial wave 
functions are expected to be roughly the same for qcqc and qqcc systems. Hence, we use the same parameter vq 
to estimate the CM interaction of + , 1 + and 2 + states of all qQqQ and qqQQ tetraquark systems with Q — c or 
b. We also use the same vq to estimate the qQQQ heavy tetraquarks. The masses are collected in Table H HTlTl As 
pointed above, the lowest-lying tetraquark state of each qQqQ system will be |QQ35, l c , 1,2) eg) \qql, l c ,0, 2) in the 
heavy quark limit, which will be too broad to be observed experimentally. Thus in Table |U we only list the masses of 
the second lowest-lying tetraquark states. The neutral ones of these states have positive C-parity. 



flavor 


(0+)V C m(vq) 


(0+)M(MeV) 


(1 + )V C m(vq) 


(l+)M(MeV) 


(2+)V C m(vq) 


(2+)M(MeV) 


qcqc 


-7.4 


3816 


-4.1 


3872 


2.7 


3988 


qese 


-6.5 


3936 


-3.8 


3982 


2.5 


4090 


scsc 


-5.6 


4056 


-3.2 


4097 


2.3 


4191 


qbqb 


-3.0 


10283 


-1.8 


10303 


0.4 


10341 


qbsb 


-2.6 


10395 


-1.6 


10412 


0.5 


10448 


sbsb 


-2.3 


10505 


-1.3 


10522 


0.5 


10553 


qcqb 


-5.4 


7046 


-3.9 


7071 


1.6 


7165 


qcsb 


-4.5 


7166 


-3.1 


7190 


1.4 


7267 


scsb 


-4.1 


7278 


-3.0 


7297 


1.4 


7371 



TABLE I: The CM energy and masses of 0+, 1 + and 2 + of qQqQ. The CM energy is in unit vq. 
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flavor 


{0+)V C m{vq) 


(0+)M(MeV) 


(1 + )V C m(vq) 


(l+)M(MeV) 


(2+)V C m(vq) 


(2+)M(MeV) 


{qq}cc 


-3.9 


3875 


1.4 


3966 


4.4 


4017 


{ss}cc 


-3.6 


4090 


0.7 


4164 


3.2 


4207 


[qq]cc 






-9.1 


3786 






qscc 


-3.8 


3982 


-7.5 


3919 


3.7 


4110 


{qq}bb 


0.7 


10346 


2.2 


10372 


3.2 


10389 


{ss}bb 


0.1 


10546 


1.4 


10568 


2.2 


10582 


[qq]bb 






-8.1 


10195 






usbb 


n a 


1 0/1/1 £ 


— 0.0 


lUoZO 


Z.D 




{qq}cb 


-1.6 


7111 


-0.2 


7135 


3.7 


7201 


{ss}cb 


-1.8 


7317 


-0.7 


7336 


2.6 


7392 


[qq]cb 


-11.1 


6948 


-9.3 


6979 


1.1 


7157 


uscb 


-9.5 


7081 


-7.8 


7110 


1.2 


7264 


TABLE II: The CM energy and 


masses of + , 1 


f and 2 + of qqQQ. The CM energy is in unit vq. [qq] is the antisymmetrical 


representation in SU(3)f, {qq} is the symmetrical representation. 








flavor 


(0 + )Vcm(vq) 


(0+)M(MeV) 


(1 + )V C m(vq) 


(l+)M(MeV) 


(2+)V C m(vq) 


(2+)M(MeV) 


qccc 


-3.5 


4979 


-4.0 


4971 


2.0 


5073 


sccc 


-3.0 


5093 


-3.1 


5091 


1.7 


5173 


qbcc 


-3.4 


8177 


-3.4 


8177 


1.3 


8257 


sbcc 


-2.8 


8292 


-2.7 


8294 


1.1 


8359 


qbbb 


-1.0 


14610 


-1.4 


14603 


0.5 


14636 


sbbb 


-0.8 


14718 


-1.1 


14713 


0.4 


14739 


qcbb 


-0.6 


11421 


-2.5 


11388 


1.1 


11450 


scbb 


-0.6 


11526 


-2.0 


11502 


0.9 


11551 


qccb 


-5.8 


8136 


-4.8 


8153 


0.9 


8250 


sccb 


-4.8 


8258 


-3.9 


8273 


0.8 


8354 


qbcb 


-4.9 


11347 


-4.6 


11352 


0.5 


11440 


sbcb 


-4.0 


11468 


-3.7 


11473 


0.4 


11543 


TABLE III: 


The CM energy and masses of + , 


1+ and 2+ of qQQQ. 


The CM energy is 


in unit vq. [qq] 


is antisymmetrical 



representation in SU(3)f, {qq} is symmetrical representation. 



IV. DISCUSSIONS 

In short summary, we have performed a schematic study of the masses of possible heavy tetraquarks using the 
color-magnetic interaction with the flavor symmetry breaking corrections. Treating X(3872) as a qcqc tetraquark and 
using its mass as input, we extract the CM interaction parameter vq for the heavy tetraquark system. With the same 
vq, we have estimated other heavy tetraquark masses. It's interesting to note from the above tables that the CM 
interaction is repulsive for 2+ heavy tetraquarks. From Table I, it's clear that + qcqc states will also exist if X(3872) 
is really a 1 + tetraquark. 
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APPENDIX A: THE RECOUPLING COEFFICIENTS 



In order to calculate the terms of V(qq), we need do the following recouplings [22, 23] 

{\qiqaD 6 (Q)D 3 (Q)S(Q)) ® \q 3 qMQ)D 3 (Q)S(Q))} [D3 . s) 
= ]T R(D 3 (Q)D 3 (Q); D 3 (13)D 3 (24); D)R(S(Q)S(Q); 5(13)5(24); S) 

x{|gifc£> 6 (13)D 3 (13)S(13)) ® \q 2 q 4 D 6 (2A)D 3 (2A)S(2A))} (D3tS) , (Al) 
where the recoupling coeffiecients are: 

R(S(Q)S(Q);S(13)S(2A);S) 

(I I S(Q)\ 

I I S(Q) , (A2) 
5(13) 5(24) 5 J 

R((^QVq)( X qVq); (Al3Ail3)(A 2 4Ai24)) 
= (-l)^+MQ+Ax3+m3 C/((10)(10)(10 )(01);(AQ M Q)(Ai3/il3)). (A3) 

We will list the recoupling coefficients of + , 1 + and 2 + heavy tetraquarks below. With qiq 2 <£> q 3 q~4 construction, 
the + heavy tetraquarks are: 

|l,l c ,5 = 0,4) = y||21,6 c ,l,2;2T,6 c ,l,2) + yi|21,3 c ,0,2;2T,3 c ,0,2) 
|405,lc,0,4> = yi|21,6 c ,l,2;2T,6 c ,l,2)-y||21,3 c ,0,2;2T,3 c ,0,2) 
|l,l c ,0,4) = y||15,3 c ,l,2;T5,3,l,2) + y||15,6 c ,0,2;T5,6 c ,0,2) 

|189,l c ,0,4) = y||15,3 c ,l,2;l5,3,l,2,)-y||15 ) 6 c ,0,2;15,6 c ,0,2) (A4) 

For + tetraquarks, the basis states of the qiq 2 (8 q 3 qi construction are: 

a° = |gi g2 21,6 c ,S'=l;M42T,6c,l);a^ = |gig221,3 c ,0;M42T,3c,0) 

a° = |qi(?2l5, 3 C , 1; g 3 g4l5, 3, 1); a° = \qiq 2 15, 6 C , 0; g3<74l5, 6 C , 0) (A5) 
The corresponding basis states of the qiq 3 <g> q 2 q4 construction are: 

b\ = \ qi q 3 l, 1,0; q 2 q A l, 1,0); b° 2 = \q iq3 35, 1, 1; «fcg 4 35, 1, 1) 

- |gi©35,8,0;M435 1 8,0);62=||«i^35 ) 8,l;( & 5435 1 8,l) (A6) 
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For the 2 + heavy tetraquarks, 



(A15) 



The transform matrix from Eq. i|A13|) to Eq. (|A14|) is: 

/ v3 _VI _\/3 V6 \ 

/ 6 6 66 3 6 \ 

Ve _\/3 Ve V3 y/e V3 

6 6 6 6 6 3 

\/3 \/t3 \/3 rj n 

3 6 3 6 u u 

\/6 V3 V6 %/3 n ri 

6 3 6 3 u " 

\/3 _\/6 _\/3 _\/6 _\/3 \/6 

6 6 6 6 3 6 

\ s/6 V3 \/6 _\/3 s/6 VS I 

\ 6 6 6 6 6 3 / 

Then the 1 + SU(6) flavor eigenstates in terms of two pairs of 9i9~ 3 <g) 9294 read 

|35,l c ,l) = i|9ig 3 l, 1,0; 529435, 1, 1) - ^|gig 3 35, 1,1; g 2 94l, 1,0) 
2 

- ol9i9335, 8,1; 929435, 8,1) — \qiq 3 35, 1, 1; 929435, 1, 1) (A16) 

3 

|280,l c ,l) = -i|9i9 3 35, 8,0; 929435, 8, l) + i|9i9 3 35, 8,1; 929435, 8,0) 

- ^ |9i 9 3 35, 8, 1; 929435, 8, l) + -|9i 9 3 35, 1,1; 9294 35, 1,1) (A17) 

6 3 

\/3 \/6 
|35,l c ,l) = — 191931, 1,0; 9294 35, 1,1) + —|9i9 3 35, 8,0; 929435, 8,1) 
3 6 

+ ^-|9i 93 35, 1,1; 9294EEO) + ^1919335, 8,1; 929435, 8,0) (A18) 
3 

/g a/3 

|35,l c ,l) = — 191931, 1,0; 929435, 1,1) —|9i9 3 35, 8,0; 929435, 8,1) 

6 3 

+ ^ |9i 9 3 35, 1, 1; 92941,1,0)- 4- 191 93 35, 8, 1;929 4 35, 8,0) (A19) 
6 3 

|35,l c ,l) = -I91931, 1,0; 9294 35, 1, 1) - -|9i9 3 35, 1,1; 92941,1,0) 

2 /2 

- -|9i9 3 35, 8,1; 929435, 8, 1) - ^-|9i9 3 3 5, 1,1; 9294 35, 1,1) (A20) 

3 6 

|280,l c ,l) = -i|9i9 3 35, 8,0; 929435, 8, l) + i|9i9 3 35, 8,1; 929435, 8,0) 

- ^ |9i 9 3 35, 8, 1;929 4 35, 8, l) + -|9i 9 3 35, 1,1; 9294 35, 1,1) (A21) 

3 



|405,l c ,2) 
= {|9i9 2 21, 6,1; 939421,6,5=1)} 

I |9i9 3 3 5, 8, 1; 929435, 8, 1) + ^19193 35, 1, 1; 929 4 35, 1, 1) (A22) 
|189,l c ,2) 

{|9i9 2 15,3, 1;939 4 15,3, 1)} 

^|9i9 3 3 5,8, 1; 929435, 8, 1) - y^|9i93 35, 1, 1; 929435, 1, 1) (A23) 

With SU C (3) and SU S (2) symmetry, similar recoupling coefficients can also be obtained. 



